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volume effects are under theoretical control for extraction of the pion charge radius. 
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I. INTRODUCTION 



Lattice QCD simulations enable first principles calculation of hadronic properties. Cur- 
rently, however, observables calculated on the lattice cannot be directly compared with 
nature due to systematic error. This error arises from approximations used in solving the 
theory numerically, such as the finite extent of the lattice, the non-vanishing lattice spacing, 
and larger than physical quark masses. For a large class of observables, effective field theories 
(EFTs) can be tailored to describe the approximations used in lattice QCD. These EFTs 
are the reliable tool to remove systematic error in extrapolating data to the physical point, 
see, e.g. With growing computational resources and improved numerical algorithms, we 
are beginning to enter a period in which lattice data in conjunction with EFT methods will 
enable first principles predictions. 

A further limitation encountered in lattice simulations is the available momentum. With 
periodic boundary conditions, the momentum modes are quantized. On currently available 
dynamical lattices, the lowest non- vanishing momentum mode is ~ 500 MeV. While EFTs 
can be tailored to address systematic error in lattice QCD, these low-energy theories cannot 
be utilized to predict the momentum dependence outside the range of their applicability. The 
extraction of quantities requiring non- zero momentum, such as phase shifts, electromagnetic 
moments and radii, is thus limited by coarse grained lattice momentum. Spatially periodic 
boundary conditions, however, are not mandated. So-called twisted boundary conditions 0, 
0, 0, 0, M, 0, 0, 0, E3, 1I| have gained recent attention because they can be utilized to produce 
continuous hadron momentum 12, 13, 14, 15, 0, 17, [HI, 19. 20l . 



The electromagnetic properties of hadrons paint an intuitive picture of the charge and 
magnetism distributions of quarks confined in strongly interacting bound states. To resolve 
such properties of hadrons, one requires momentum transfer and this presents a challenge 
for the lattice. Calculations of the pion form factor from lattice QCD, for example, have 
matured since the original pioneering works 21, 22|, but still suffer from systematic error. 
Recently it was demonstrated [23[ that current dynamical calculations (24 . I25L 26 ] of the 
pion charge radii are limited by systematic error relating to the momentum, chiral, and 
possibly also continuum extrapolations. Here we show that the uncertainty surrounding 
the momentum extrapolation can be eliminated with flavor twisted boundary conditions. 
With these imposed, one can extract the charge radius from simulations at zero Fourier 
momentum. 

As boundary conditions modify the l ong -range physics on the lattice, we must be careful to 
address effects of the finite volume. In |14j volume effects for single-particle matrix elements 
with twisted boundary conditions were shown to be similar to those with periodic boundary 
conditions. Thus, provided one remains below multi-particle cuts, the effect of the finite 
volume is exponentially suppressed in asymptotic volumes. In practice, lattice volumes 
are not asymptotically large, and the study with partially twisted boundary conditions 
in 20| found sizable volume corrections to the nucleon isovector magnetic moment, ~ 25%, 
even for m^L « 4. The bulk of such corrections arises from finite volume terms that 
do not respect hypercubic invariance. Physically one expects contributions of this type 
because twisted boundary conditions introduce a particular direction in space. In this work, 
we determine hypercubic breaking corrections to the pion dispersion relation for partially 
twisted boundary conditions. These corrections result in a finite volume modification to 
the pion's momentum. Theoretical control over these volume effects is necessary wherever 
twisted boundary conditions are employed because of the dominance of long-range pion 
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physics in low-energy hadronic observables. Additionally we determine the finite volume 
effects for the extraction of the pion charge radius from lattice simulations with partially 
twisted boundary conditions. We further demonstrate numerically that such corrections are 
on the order of a few percent on current dynamical lattices. 

Our presentation has the following organization. In Section [TTJ we briefly review partially 
twisted boundary conditions and their incorporation into partially quenched chiral pertur- 
bation theory (PQ%PT). By calculating the pion self-energy, we determine the form of the 
pion dispersion relation in a finite box and demonstrate that the induced momentum due to 
twisted boundary conditions is modified by volume effects. The isopsin splitting among the 
pions is also evaluated in PQxPT. Next in Section ITTlj we show that the pion electromagnetic 
form factor can be determined on the lattice at continuous values of the momentum trans- 
fer. Using PQxPT, we then determine the finite volume corrections resulting from partially 
twisted boundary conditions. These volume effects are a controlled theoretical error that 
must be removed in order to isolate the pion charge radius from lattice data. A summary 
(Section [TVT) concludes the paper. 

II. DISPERSION RELATION WITH TWISTED BOUNDARY CONDITIONS 

The quark part of the partially twisted, partially quenched QCD Lagrangian reads 

£ = Q($+m Q )Q, (1) 



where the six quark fields have been accommodated in the vector Q = (u,d,j,l,u,d) T 
that transforms in the fundamental representation of the graded group S77(4|2). The 
mass matrix in the isospin limit of the valence and sea sectors is given by uiq = 
diag(m u ,m u ,mj,mj,m u ,m u ). The valence u and d quarks are accompanied by equal mass 
bosonic ghost quarks u and d. Functional integration over both valence and ghost quarks 
thus produces no net determinant factor, while integration over the sea quarks j and / pro- 
duces a determinant. The graded symmetry of SU(4\2) hence provides a way to write down 
a theory corresponding to partially quenched QCD. For a thorough discussion, see [321 ]. 

All Q fields appearing in the Lagrangian satisfy periodic boundary conditions, and the 
effect of partial twisting has the form of a uniform U(l) gau ge field: = + iB^, with 



the flavor matrix B^ = diag(-B^, B®, 0, 0, B™, B^), see [14|, [16| for details. The field B^ acts 
as a flavor dependent field momentum: for the u quarks B™ = (0, 6 U /L), and similarly for d. 
The twist angles can be varied continuously and cannot be renormalized by short distance 
physics. Notice that the twists are not implemented in the sea sector of the theory; thus, 
new gauge configurations do not need to be generated for each value of 6. 

The dynamical effects due to partial twisting can be addressed using the low-energy ef- 
fective theory of partially quenched QCD in a finite volume. This low-energy theory is 
PQxPT which describes the dynamics of the Goldstone bosons that arise from spontaneous 
breaking of chiral symmetry. We determine volume effects in the p-regime of chiral per- 
turbation theory [28| in which the longest-range fluctuations of the Goldstone modes do 
not restore chiral symmetry in a finite volume. The effective theory PQxPT is written in 
terms of a periodic coset field £ = exp(2i(f>/ f) into which the meson fields 4> are embedded 
non-linearly. To incorporate partially twisted boundary conditions, we couple the field S to 
the background gauge field B^. In terms of the field E, the Lagrangian of PQxPT appears 
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FIG. 1: Pion self-energy diagrams in PQxPT. Dashed lines represent mesons, while the cross 
represents the partially quenched hairpin interaction. 



as [23, |29|, [3d, |31|, |32j 



C = ^str [b^D^) - A str (m^E + E^q) + (2) 

With /i — »■ oo, the flavor singlet field $ — str(E)/v2 is integrated out, but the resulting 
flavor neutral meson propagators have both flavor connected and disconnected (hairpin) 
contributions 32|. The action of the covariant derivative t)^ is specified by 

= c^E + i[B„, E] + %Al [T a , E] , (3) 

where we have additionally included an isovector source field which will be utilized 
below in Section HTT1 The block diagonal form of the super-isospin matrices T a is chosen to 
be T a = (T a ,0,T a ), where T a are the usual isospin generators 20|]. From the commutator 
structure that couples the background gauge field to E, one sees that mesons have a 



flavor- dependent momentum boost [14j. Expanding £ to tree level, one finds that mesons 
with quark content QQ' have masses 1 

m QQ' = T^( m Q + m Q')> ( 4 ) 

and kinematic momenta 

Pqq> =P + B q -B q/ , (5) 

where p is the Fourier three-momentum. 

Beyond tree-level, observables receive non-analytic corrections from loop diagrams gen- 
erated in the effective theory. These loop diagrams, moreover, are modified in finite volume. 
As twisted boundary conditions are flavor dependent, isospin symmetry is broken and the 
degeneracy among the pions is lifted by volume effects. The isospin splitting among the pi- 
ons was determined for completely twisted in [13]. Here we determine the splitting in 
partially twisted PQxPT. Moreover as twisted boundary conditions introduce a particular 
direction in space, the pion dispersion relation is no longer a hypercubic invariant function 
of the Fourier momentum p. We determine the form of the dispersion relation and show 
that the kinematic momenta in Eq. §5§ are consequently modified. 

To determine these volume effects, we must calculate the one-loop diagrams contributing 
to the self-energy of the pions. These diagrams are depicted in Fig. [lj Evaluating the 
contractions leading to the diagrams in the Figure, we are confronted with various terms 
that can be grouped according to where the two derivatives of the four-pion vertex act. Let 
us detail the calculation of the diagrams for the charged pion. Terms where both derivatives 



1 We shall often write m^, and m^o to denote the mesons formed from purely valence quarks. Other cases 
for which we also use as an abbreviation will be clarified in the text. 
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at the vertex act on the internal line or both act on the external legs are the most familiar 
because these survive the infinite volume limit. The former give rise to a mass shift, and 
the latter constitute a wavefunction renormalization correction. A general contribution to a 
loop diagram in finite volume, Cl, can be written as 

C L = C O0 + 8 L {C), (6) 

where is the result for contribution C in infinite volume and the difference Sl(C) = 
Cl — Coo is ultraviolet finite. Combining the mass shift and contributions from wavefunction 
renormalization, we arrive at <5z,(m^±) the finite volume shift of the charged pion mass which 
is given by 

h ^ =$Jk [(ml " m - )Jl/2( °' ^ • (7) 
As is well known, due to delicate cancellations the charged pion mass only receives loop 
contributions from flavor disconnected (hairpin) propagators. Thus so too must the finite 
volume correction. These hairpin contributions are most economically written in terms of 
derivatives with respect to the pion mass squared, as appears in Eq. ([7]). Because the net 
contribution to the mass is from flavor neutral loop mesons, there is no -B-dependence in 
the volume effect at this order. The function X ( g(S,m 2 ) results from taking the difference 
of the finite volume contribution minus the infinite volume contribution, and is given by 

. 2 1 x ^ 1 r dk 1 

f>{ > m ) = 7/ 2^ + B y + m 2y - J (27i) 3 [(k + B) 2 + mY' ^ 

To evaluate this finite volume function, the sum over Fourier momentum modes can be cast 



into an exponentially convergent form by using the Poisson re-summation formula. Resulting 



expressions are then swiftly evaluated using Jacobi elliptic functions, see [14], |23j, [33 . 

We are not done evaluating the momentum contractions for the diagrams shown in the 
Figure. There are also contributions where one derivative acts on the internal line and one 
derivative acts on either of the external legs. These contractions yield terms proportional to 

1 y (p + B„±) • (k + B q ) 

^ 3 fc J(k + Bj^ q ' 

where q is a valence quark index and m| g is the mass of the corresponding valence-sea meson 
formed from q and either of the degenerate sea quarks. As such contractions involve only 
valence-sea meson loops, they are absent the quenched approximation. Notice that in the 
infinite volume limit, the sum in Eq. (Q can be converted into an integral. Using a shift of 
variables, the infinite volume contribution is proportional to 

dk (p + B w ±) ■ k 



( 27r ) 3 + m% 



and vanishes by 5*0(3) invariance. Furthermore in a lattice with periodic boundary condi- 
tions (B u = Bd = 0) the sum in Eq. (jSJ) becomes 
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FIG. 2: Maximal change in the induced momentum due to twisting. Here, the pion mass m n is 
the valence-sea meson mass, and U = — 6 d in the limit of vanishing twist. 

and vanishes by cubic invariance. With twisted boundary conditions, however, the sum 
in Eq. ([9]) does not vanish. Hence terms that ordinarily vanish in infinite volume due 
to rotational invariance (and similarly vanish in a periodic box due to cubic invariance) 
contribute. The sum, moreover, is non-vanishing only in the directions for which B q is 
non-vanishing. These contributions then are clearly due to the breaking of cubic invariance 
introduced by the spatial direction specified by the background gauge field. To describe 
finite volume effects from these terms, we define the function 

which can be evaluated in terms of a derivative of the elliptic Jacobi function [2p| . 

Combining all effects due to the finite volume, we find that the charged pion dispersion 
relation has the form 2 

E^(p? = (p + B n± ± Kf + ml + 6 L (ml ± ) , (11) 

where the function K appears as an additive renormalization to the field momentum B n ± = 
±(B U — Bd), and is given by 

K = -l[K 1/2 (B U7 m%) - K 1/2 {B d ,m%)]. (12) 

From Eq. (jl2]) . we see that the pion field momentum is renormalized only in the direc- 
tions for which B n ± is non-vanishing. Furthermore, the analogous quenched finite volume 
contribution vanishes. 

We can investigate the finite volume renormalization of the field momentum by consid- 
ering the relative change in the pion's momentum due to twisting 

AB„ + /B w+ = J2 Kj / B i+- ( 13 ) 

3 



2 In order to arrive at Eq. (JTTJ) , we added a higher-order term proportional to K 2 which scales as / 4 . The 
one-loop results scale as f~ 2 . 
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FIG. 3: Relative change in the induced momentum due to twisting as a function of the twist angle. 
The pion mass m n refers to the valence-sea meson mass, and the box size is fixed at 2.5 fm. 

The relative change in the negatively charged pion's momentum is numerically identical. In 
the limit B n + —>■ 0, we have AB n + /B 7T + — Vb ■ K. One should keep in mind, however, that 
neither of these quantities is accessible without twisted boundary conditions, but the twisted 
divergence serves as an approximate relative change for small twists. In Fig. [21 we plot the 
relative change in induced momentum AB n + / B n + in Eq. f[T3"j) for the worst case scenario. We 
take 6 U = —6 d and further choose each component near zero. Finite volume modifications for 
this case are seen to be substantial from the figure. Thus in applications of twisted boundary 
conditions, one must take care to avoid large additive volume renormalization of the twisting 
parameters. Choosing only one flavor to be twisted, e.g. d = 0, and implementing the twist 
in only one spatial direction reduces the volume effect by a factor of six. Furthermore in the 
actual implementation of twisted boundary conditions near vanishing twist angles will not 
be employed. Thus in Fig. [31 we investigate the behavior with respect to the single twist 
angle 9, where U = (0, 9, 0) and 6 d = 0. The relative change in the induced momentum 
is plotted versus 9 for a few valence-sea pion masses (denoted by m n in the figure), but 
with the box size L = 2.5 fm fixed. The behavior with 9 is damped oscillatory and crosses 
zero when 9 = (2n + l)ir. For reasonably small values of the twist angle, 9 ~ 7r/8, and at 
(valence-sea) pion masses ~ 0.3 GeV, the induced momentum is shifted by ~ 2% from volume 
effects. Thus the finite volume renormalization of the twist angles can practically be ignored 
on current lattices. As pion masses are brought down, however, volume renormalization of 
the induced momentum may need to be accounted for by evaluating Eq. f[T2l . 

The above discussion applies to the charged pions 7T ± . The neutral pion dispersion rela- 
tion, however, is unaffected by twisting. Momentum contractions where one derivative acts 
on the internal line and one acts on the final state are canceled by those where the second 
derivative acts on the initial state. Thus only the mass of the neutral pion is affected by the 
finite volume. Performing the self-energy calculation for the 7r°, we find the finite volume 
mass shift 




Combining the finite volume shifts for the charged and neutral pions, we find that the 
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FIG. 4: Maximal isospin splitting between the charged and neutral pions. 

isospin splitting due to partially twisted boundary conditions has the form 

Am 2 = 5 L (ml±) - 5 L (m 2 ) 

9 

= ^f[l 1/2 (0,ml)~l 1/2 (B n+ ,ml)}, (15) 

where the hairpin contributions have canceled out. The splitting Am 2 is only sensitive to 
the valence pion mass and vanishes when the twists preserve isospin, B u = B&. In Fig. HI we 
plot the relative isospin splitting Am 2 /m 2 as a function of the lattice size L for a few values 
of the valence pion mass m n . The maximal isospin splitting occurs for 6^+ = (tt, tt, it), and 
we have chosen to plot this value to investigate the worst case scenario. From the figure, 
we see ~ 10% maximal isospin splitting at the physical pion mass in a relatively small box. 
For lattice pion masses ~ 0.3 GeV, the maximal isospin splitting is ~ 2% in a 2.5 fm box. 
For applications of twisted boundary conditions, such as the extraction of the pion charge 
radius, only one component of 6 n + need be non-vanishing. Moreover, this twist angle should 
be chosen less than it to reach smaller momentum transfer. These considerations restrict 
the isospin splitting to be < 1% for pion masses ~ 0.3 GeV in a 2.5 fm box. For this reason 
we neglect the isospin splitting in virtual loop meson masses below. 

Having detailed the finite volume corrections to the pion mass and dispersion relation, we 
now describe how partially twisted boundary conditions enable the extraction of the pion 
electromagnetic form factor. 



III. PION ELECTROMAGNETIC FORM FACTOR 



The pion electromagnetic form factor G n (q 2 ) arises in the charged pion matrix element 
of the electromagnetic current J^ m , namely 



(n + (p')\J™\n+(p)) = (p' + p) fl G 7T ( q 2 



(16) 



with q = p' — p. In the isospin limit, however, vector flavor symmetry relates this matrix 
element to the isospin transition 20[ 

<7T + ( P ')i^kV)> = -v2> + (p')urk>)>> ( i7 ) 



s 




FIG. 5: Contraction for the pion form factor in Eq. (|17p . The source and sink are denoted by filled 
circles, the operator insertion by an open circle. Single lines represent twisted d-quark propagators, 
while the double line represents the twisted ti-quark propagator. 



where we have made use of charge conjugation invariance, and the phase convention is 
that commonly employed in xPT. Thus by calculating the pion form factor from the left- 
hand side of the above equation, we can induce momentum transfer with partially twisted 
boundary conditions The neutral pion is flavor non-singlet and because we work in 

the isospin limit, the annihilation diagrams contributing to the left-hand side of Eq. (TP?]) 
exactly cancel. Thus only the quark contraction shown in Fig. [5] must be calculated in 
lattice simulations. While flavor neutral states in partially twisted (and partially quenched) 
theories have sicknesses related to lack of unitarity the illness does not infect the 



states with isopsin 1 = 1. Consequently the propagator of the it is diagonal with a simple 
pole form, and asymptotic tt° states can be created a la LSZ. 

Application of partially twisted boundary conditions to the matrix element in Eq. ([1] 
thus leads to the kinematical effect 

(7r + (0)|n 74 d|7r°(0)> = -V2(p 4 + p,)G^B 2 ), (18) 

where for emphasis we have projected the initial and final states onto zero Fourier momen- 
tum, and restricted our attention to the fourth component of the current. The energy of 
the initial state is p 4 = im n and the energy of the final state is p' A = iy/m% + B 2 + , while 
the momentum transfer squared is B 2 = 2m 2 (— 1 + yjl + B 2 + /m^j . Finally the induced 
three-momentum transfer B^+ is given by the difference of twists of the flavors changed 
B n + = B u — Bd- In writing the above expression, we have temporarily ignored dynam- 
ical effects from the finite volume. To assess the effects of the finite volume on the pion 
electromagnetic form factor with partially twisted boundary conditions, we must again use 
PQxPT and do so in the p-regime. 

At next-to-leading order in the EFT, there are local contributions from higher-order 
operators and loop contributions generated from the leading-order vertices. We handle the 
local contributions first. In partially twisted PQ^PT, the isovector current operator J+ 
receives a local contribution at next-to-leading order from the term 

5J+ = ia 9 D v str [t+ D U E] - [D V E\ D^E] 

This is just a rewriting of the current derived from the «g term of the Gasser-Leutwyler 



Lagrangian [34]. As explained in 20], the derivatives appearing in operators for external 
currents must be upgraded to .D's. There are also local contributions to the current from 
the CK4 and a*, terms of the fourth-order Lagrangian, but these are canceled by wavefunction 



renormalization, see 35, 36 



The one-loop contributions to the isovector form factor are shown in Fig. [6] (along with 
the requisite wavefunction corrections shown in Fig. [1]) and lead to finite volume modifica- 
tions. The hairpin diagrams exactly cancel in both infinite and finite volumes. This can 
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FIG. 6: One-loop contributions to the isospin transition matrix element in PQxPT. The wiggly 
line shows the insertion of the isospin raising operator. Other diagram elements are the same as 
in Fig. [TJ 



be understood easily due to the isospin relation in Eq. (TT7l) . On the right-hand side of this 
relation is the matrix element of the electromagnetic current. Contributions to this matrix 
element at one-loop are only due to charged meson loops, thus there can be no hairpins which 
arise from flavor (and thus charge) neutral mesons. Explicit calculation of the left-hand side 
of Eq. (fTTj) verifies this. 

Returning to Eq. (jTSI) but now including the effect of the finite volume at one-loop order 
in the chiral expansion, we have 



(tt+CpOI U 7 4rf|vr°(p)) = -V2i^[E w+ (p') + EAP)] [G^{Q 2 ) + G fv ] 

+ [E n+ (p') - EAP)] [Gf v + (p + B w+ +p) ■ Gf v ] j, (19) 

where E^+ijp') is given in Eq. (ITT]) , and E^o^p) 2 = p 2 + m 2 + (^(m 2 ^). Above Q is the 
effective momentum transfer, Q = q + B^, and G^jQ 2 ) is the partially quenched pion 
electromagnetic form factor in the infinite volume limit [37], namely 



GM 



409^2 1 

f 2 ^ 48vr 2 / 2 



Q 2 log^ + 4m^ 



-Q 2 



where the auxiliary function jF(x) is defined by 



(20) 



x 



-log 

X 



1 



- + ie 

X 



= + -x-2. 

1 - 1 + is + 1 6 



(21) 



The remaining contributions in Eq. (TT9~|) are finite volume modifications and have been 
denoted so with a FV subscript. There are various flavor and momentum contractions for 
the meson fields that lead to these terms. They fall into three categories. We give the 
explicit results for each, and comment on the physics behind such terms. 
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The first finite volume contribution Gfv is the finite volume modification to the form 
factor. This contribution arises from the difference between the sum of Fourier modes in the 
loop and the infinite volume loop integral that leads to Eq. (120]) . This effect remains even 
in the absence of twisted boundary conditions. Explicitly we have 



G 



FV 



2P 



dx 



Z1/2 {xQ - B u , m 2 u + x(l - x)Q 2 ) - li /2 (-B u , m 2 u ) 



+ J1/2 (xQ + B d , m 2 ju + x{\ 



x)Q 2 ) -X 1/2 (B d ,m%) 



(22) 



and when B u = B d = 0, we recover the finite volume modification to the form factor with 



periodic boundary conditions 23 



The remaining finite volume effects in Eq. (fT9|) are due to symmetry breaking introduced 
by the partially twisted boundary conditions. The second finite volume contribution G F y 
arises from the breaking of isospin and vanishes when B u = Bd- This term moreover enters 
the matrix element in Eq. (jT9j) proportional to the energy difference between the initial and 
final states. Ordinarily such an additional structure is forbidden by current conservation. In 
finite volume, however, the partially twisted boundary conditions break isospin symmetry 
and the isospin current J+ is no longer conserved. This finite volume correction G F ° V appears 
as 



6/ 



l — ^ dx{l -2x)|4J 1/2 (xQ,ml + x{l - x)Q 2 ) 

-Q 2 [1 + 2x(l - x)} X 3/2 {xQ, m 2 n + x(l - x)Q 2 ) 



1 



+ (1 - 2x)Q ■ 7C 3/2 (xQ, ml + x(l - x)Q 2 ) \ - -Am 2 /m 2 , (23) 

and involves the isospin splitting among the pions Am 2 given in Eq. ffTBT) . 

The final finite volume modification G v p V arises from the breaking of hypercubic invari- 
ance introduced by partially twisted boundary conditions. This contribution is given by 



Grot 



FV 



i r 



IC3/2 {xQ - B u , m 2 u + x(l - x)Q 2 ) 



+ /C3/2 (xQ + B d ,m 2 ju + x(l 



x)Q 2 ) 



(24) 



Notice this term only vanishes when B u = B d = 0. When B u = B d 7^ 0, isospin symmetry 
is restored and we might expect G r p V to vanish given the way this term enters Eq. f|T9l) . 
This is not the case, however, because twisted boundary conditions break cubic invariance 
and extra terms are allowed. Indeed the decomposition of the current matrix element in the 
case when the twisted boundary conditions preserve isospin B u = B d 7^ 0, 



(n+(pi)\J+\n\p)) = (p>+p),G(Q 2 ), 



(25) 



relies on 5*0(4) rotational invariance. While this reduces to hypercubic invariance on the 
lattice, if there are no terms in the EFT that violate 5*0(4) at this order, then Eq. fl25|) 
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still holds. Twisted boundary condition, however, do not respect the hypercubic invariance 
of the lattice and the above form is no longer valid. Explicitly G T p V is not invariant un- 
der the transformation Qj — > —Qj, and indeed breaks hypercubic invariance even in the 
isospin symmetric case. We should thus interpret this term in Eq. (1191) as arising from cubic 
symmetry breaking rather than the violation of isospin. 3 

Having deduced the volume corrections to the pion form factor using PQ^PT, we can 
now investigate the impact of partially twisted boundary conditions on the extraction of 
the pion charge radius. To do this, we choose for simplicity B d = and B u = (0, |B|,0), 
with \B\ = 8/L. Furthermore we take the source and sink to be projected onto zero Fourier 
momentum. We investigate the relative difference of the form factor (minus the isovector 
charge) in finite volume to infinite volume. Specifically we investigate the following quantity 

1 (tt+(o)|j 4 + | 7r°(o)) - GAQ 2 ) 

S L [G W (Q 2 ) - 1] = ^ (P4+P4) GAQ2) _ i . (26) 

For small twist angles, <C 1, the relative difference in the form factor 5l[G. k (Q 2 ) — 1] 
becomes the relative difference in the charge radius 

limS L [G w (Q 2 )-l]=S L [< rl >} 

a — >\J 



<rj> L - <r| >oq 



(27) 



Finally as expressions depend on both the valence-valence and valence-sea meson masses, 
we take the theory to be unquenched, so that m 2 u = m 2 . In Fig. [7J we plot Sl[< r 2 >] 
as a function of L for various values of the pion mass to contrast our results with those of 
Ref. [23| . We plot the total contribution which arises from the three finite volume terms in 
Eq. (fTPl . however, nearly all of the effect is due to Gpv- Contribution from terms that arise 
from the breaking of symmetries in finite volume, G F ° V and G r py, are suppressed relative 
to Gpv due to the kinematical pre-factor accompanying these terms. The mag nitude of 
the finite volume effect shown in the figure is smaller than that calculated in [23(, and has 
differing sign. In that study, however, periodic boundary conditions were assumed and the 
extraction of the charge radius required utilizing the smallest available Fourier momentum 
transfer. With twisted boundary conditions, we have eliminated this restriction. Volume 
effects are consequently different. 

Of course to extract the charge radius from twisted boundary conditions, we require 
6 7^ 0, and thus we need to investigate SlIG^^Q 2 ) — 1] in Eq. (12"5"|) . In Fig. [HI we fix the 
lattice volume at 2.5 fm and plot the relative difference SlIG^IQ 2 ) — 1] as a function of 6 
for a few values of the pion mass. We see that the finite volume effect increases with 9 
over the range plotted. Generally the finite volume effects for form factors decrease with 



momentum transfer [201 . [38|, because as the resolving power of the virtual probe increases, 
the effect of the finite volume decreases. The behavior with increasing momentum transfer 
(here parametrized by 6) is indeed damped oscillatory but the oscillations appear beyond 



As for the other finite volume contributions, Gfv m Eq. (|22l) is not invariant under Qj — > —Qj unless 
B u — Bd- The contribution G^y is invariant under Qj — > — Qj and so its presence in Eq. (fl9|) is 
unambiguously due to isospin violation. 
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FIG. 7: Finite volume shift of the pion charge radius with partially twisted boundary conditions. 
Plotted versus the box size L is the relative change in the charge radius 5l[< r% >] for a few values 
of the pion mass. 

the range of the EFT. Nonetheless, the finite volume effect only amounts to a few percent 
on current dynamical lattices. 

Flavor twisted boundary conditions introduce systematic error into the determination of 
observables, such as the pion charge radius. We stress, however, that this is a controlled 
effect. There are no extraneous phenomeno logical form factor fitting functions needed. The 
EFT can be used for the momentum extrapolation because employing 9 < n brings us into 
the range of applicability of the EFT on current dynamical lattices. Furthermore there are 
no additional parameters introduced in the EFT to describe the volume dependence with 
twisted boundary conditions. For mesonic observables, knowledge of the pion decay constant 
fn and the pion mass m n (possibly also the valence-sea pion mass) is enough to remove the 
volume dependence. Higher-order dependence on the quark mass, volume and momentum 
transfer can all be calculated systematically in the EFT. In particular, the pion charge radius 
determined from twisted boundary conditions will be rather insensitive to volume effects on 
current dynamical lattices, see Figs. [7] and [SJ 

IV. SUMMARY 

Above we have shown that in the isospin limit the pion electromagnetic form factor may 
be extracted from the isospin changing matrix element in Eq. ( TT7|) . Utilizing flavor twisted 
boundary conditions allows one to circumvent the restriction to quantized momentum trans- 
fer. We have detailed the finite volume corrections to the pion mass and dispersion relation 
using partially twisted PQxPT. Addressing these effects is crucial to determining observ- 
ables, such as the pion charge radius, in simulations with twisted boundary conditions. 

We find that the twist angles acquire additive finite volume renormalization which can 
be sizable near vanishing twist. This effect is only present for dynamical simulations, and 
the renormalized twist angle 6 r q en (L), with the quark flavor q = u or d, is given by 




(28) 
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FIG. 8: Finite volume shift of the pion form factor with partially twisted boundary conditions. 
Plotted versus 9 is the relative change in the form factor 6l[G 7T (Q 2 ) — 1] for a few values of the 
pion mass. The three-momentum transfer from twisting is \B\ = 9 x 0.079 GeV, and the box size is 
fixed at 2.5 fm. 

where j labels the sea quark (in a two- flavor degenerate sea). For practical applications of 
partially twisted boundary conditions on current lattices this volume correction is ~ 2%. In 
finite volume, partially twisted boundary conditions lead to isospin symmetry breaking. The 
volume induced isospin splitting among the pions is under good theoretical control, < 1% 
for practical applications of twisted boundary conditions on current lattices. Finally, the 
extraction of the pion charge radius with partially twisted boundary conditions introduces 
no new free parameters, and finite volume effects can be assessed with effective field theory 
methods. We find that the finite volume effects for extracting the pion charge radius from 
simulations with twisted boundary conditions are only a few percent. 

A simple exercise in isopsin algebra j22|, 23[ shows that the form factors of all isospin 



1 = 1 mesons (pseudoscalar, vector, . . .) can be computed on the lattice at continuous 
values of the momentum transfer. Thus, for example, the magnetic moment and electric 
quadrupole moment of the p + meson are accessible at zero Fourier momentum. Apart from 
the pseudoscalar mesons, the volume effects cannot be calculated systematically, however, 
as there are no known EFTs for these mesons with a controlled power counting. One 
would have to demonstrate that the electromagnetic properties extracted from the lattice 
are stable under changes in the volume. Nonetheless simulating lattice QCD with flavor 
twisted boundary conditions provides a promising technique to overcome the restriction to 
quantized momentum transfer. 
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